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SUMMARY 


Graphs of phase and group refractive indices computed from the collision- 
less magnetoionic theory are presented. The indices are computed for both 
the entire range of electron concentration and the intensity of the earth's 
magnetic field to be encountered by operational topside ionospheric sounders 
in the International Satellites for Ionospheric Studies (ISIS) series 
(Alouette I, Explorer XX, Alouette II, ISIS-A, ISIS-B) . 


INTRODUCTION 


Swept high-frequency radar soundings of the earth's ionospheric plasma 
have been performed routinely and bottoms ide ionograms collected for over 
40 years by more than 100 ground-based stations. Topside ionograms obtained 
from sounders on board the earth-orbiting satellites Alouette I and 
Alouette II have been collected since 1962 . The theory of propagation of 
electromagnetic waves in an ionized medium (i.e., the magnetoionic theory) is 
well documented and generally accepted (refs. 1 and 2). The key parameters 
affecting the propagating waves are the phase and group refractive indices 
which, in turn, define the wave phase and group velocities relative to the 
velocity of light in free space. These indices are functions of the radar 
wave frequency, the local ionospheric electron concentration, the earth's mag- 
netic field intensity, the angle between the directions of the wave normal and 
magnetic field, and the particle collision frequencies. Computations of 
phase and group refractive indices appropriate to the ionosphere below the 
F 2 layer peak, for restricted ionospheric models and for particular geograph- 
ical locations, have been published (refs. 3 - 5 )* The earth-orbiting radar 
sounder has made it necessary to consider these indices at substantially all 
geographic longitudes and latitudes at heights from 300 to 3000 km. Within 
this range the collisionless magnetoionic theory is applicable. It is the 
purpose of this report to present graphs of the indices appropriate to the 
entire range of electron concentration observable during a sunspot cycle 
throughout the above geographical region. Applications of these results to 
the reduction of topside sounder data are reviewed by Jackson (ref. 6 ). The 



results are also applicable to other types of experiments involving propaga- 
tion of electromagnetic energy through cold magnetoionic plasmas. 


REFRACTIVE INDICES 


The phase and group refractive indices are defined by 


h = c/v p 

(1) 

■ ' = c/v g = p + f(dp/df) 

(2) 


where 


p phase refractive index 

p' group refractive index 

v phase velocity 

Jr 

v group velocity 

S 

c velocity of light in free space 

f electromagnetic wave frequency 

The collisionless phase refractive index is given by the 
equation of the magnetoionic theory (ref. l): 


Appleton-Hartree 
11/2 


[X = 


1 


X 


1 _ jj + /ll_sinl_e + Y 2 cos 2 q 

2(1 - X) Jb(l - X) 2 


( 3 ) 


where 

X = f u 2 /f 2 = N/l2 , k-OOf 2 
N electron concentration, electrons/cm 3 

f electromagnetic wave frequency, Mc/s 

f plasma frequency, Mc/s 

and 

Y = f h /f = 2 .8 b/ f 
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B 

% 
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magnetic induction , gauss 
electron gyrofrequency, Mc/s 

angle between the directions of the wave normal and earth's magnetic 
field 


It can be seen that two modes of propagation are possible in a magnetoionic 
medium, the ordinary mode (with the plus sign) and the extraordinary mode 
(with the minus sign) • These are usually denoted with appropriate subscripts 
x and o. For clarity the ordinary and extraordinary phase refractive indices 
may be written: 

\l/2 
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Ho = 1 - 


rj 


where 


hx 




1/2 


S„ = 1 - 


Y 2 sin 2 0 YR 

+ 


2(1 - X) 2(1 - X) 


(5) 


(6) 


S-y — 1 


Y 2 sin 2 0 


YR 


2(1 - X) 2(1 - X) 


(7) 


R = [Y 2 sin 4 0+4 cos 2 0(1 - X) 2 ] 


1/2 


( 8 ) 


Applying equation (2) to equations (4) through (8) yields the ordinary 
and extraordinary group refractive indices: 


Ho 



XY(1 - X)cos 2 0 r _ Y sin 2 0 / j + xY 
S 0 2 (Y sin 2 0 + R) L R V 1 “ V. 


4 = - 

X H x 


1 + 


X(1 


S x ) 


2S X 2 


1 + 


Y sin 2 0 fl + X 
1 - 


(9) 


( 10 ) 


These equations simplify considerably for the special cases of longitudinal 
propagation (0 =0°) and transverse propagation (0 = 90° ): 
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(a) longitudinal propagation 


Ho = 1 - 


Hx = 1 ~ 


1 + Y 


1 - Y. 


1/2 


1/2 




1 - 


XY 


2(1 + Y) 2 J 


^ = - 
x I 

(b ) transverse propagation 


1 + 


XY 


^ x H x r ’ 2(1 - Y) 2 J 


Mo = (1 - X) 


1/2 


Mx 


1 - 


X(1 - X) 

1 - X - Y 2 J 


1/2 


< - r 
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= 1 


1 + 


XY 2 


x M X L (1 - X - Y 2 )* 
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( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


GRAPHS OF p AND p' 


The variation of p 2 with X for a given value of Y is compactly- 
illustrated in figure l(a) for Y < 1 and in figure l(b) for Y > 1. The 
vertical cross-hatching bounded by the solid curves represents the computa- 
tion of p Q 2 from equation (4) while the horizontal cross-hatching bounded 
by the dotted curves represents the computation of p x 2 from equation (5) . 
The longitudinal and transverse limits are labelled L and T, respectively, 
and are computed from equations (ll) , (12), (15) > and (l6) . The specific 
case of Y = l/2, 0 = 45° is plotted in figure l(a) and the case of Y = 3> 

0 = 15° is plotted in figure 1(b). 

For Y < 1 (fig. 1(a)) there are one branch for the ordinary mode 
(0 ^ X £ l) and two branches for the extraordinary mode (0 ^ X ^ 1 - Y and 
(1 - Y 2 )/ (l - Y 2 cos 2 0) £ X £ 1 + Y) . The latter branch of the extra- 
ordinary mode is referred to as the Z mode. For Y > 1 (fig* 1(b)) the 
ordinary mode has two branches (0 £ X £ 1 and (l - Y 2 )/(l - Y 2 cos 2 0) £ X; 
the latter branch is called the whistler mode and exists only when 
Y 2 cos 2 0 >l) while the extraordinary mode has only one branch, the Z mode 
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branch (0 ^ X ^ 1 + Y). The ordinary , extraordinary, and Z modes are 
observed as distinct reflection traces in topside ionograms . 

Zeros occur in p at X = 1 (ordinary mode), X = 1 - Y (extraordinary 
mode), and X = 1 + Y (Z mode). An infinity occurs in p at 

X = (l - Y 2 )/[l - Y 2 cos 2 0] (z mode for Y < 1 and ordinary mode for Y >l). 

Both the zeros and infinities in p result in infinities in p* . It should 
be noted that the infinity in p' associated with X = (l-Y 2 )/[l-Y 2 cos 2 0] 

is not a wave reflection condition as are the other infinities . 

For the extraordinary mode when Y > 1 (i.e., the Z mode), p and p* 
exhibit discontinuities at X = 1 for 0=0 (longitudinal propagation) , and 
p* exhibits a maximum at X = 1 for other angles. It can be shown (ref. 3) 
that the following relationships obtain: 



v \ l/2 
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The variations of p' vs. X for different values of Y and 0, plotted 
on a reciprocal scale to exhibit the infinities, are shown in figures 2 
through 113. Table I contains a key to the order of these figures. Figures 2 
through 66 are plots of p' vs. X for fixed 0, variable Y. Figures 67 
through 113 are plots of p* vs. X for fixed Y, variable 0. 


TABLES OF p AND p'. 


For the reader who requires greater accuracy than is obtainable from 
figures 1 through 113 > tabulations of p and p* may be obtained by mailing 
the request card in the back of this report. The tabulated quantities 
include P Q 2 , p x 2 , Pq, p x , S x and R. The range and increments of the 

independent variables are: 


X 

0(0. 1)1. 9 + Y 

for 

Y s 1 


0(0.1)10 

for 

Y > 1 

Y 

0(0. 1)5 



0 

o°(io 0 )90° 

and 

15°; ' 
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When values of X or Y other than these are encountered, it is suffi- 
cient, for most purposes, to interpolate between the tabulated results. 
Interpolation for values of Q other than those given may not be sufficiently 
accurate, however. A Fortran IV listing of a program to compute the phase 
and group refractive indices for any Q and for the X and Y intervals and 
increments mentioned above is given in table II. 
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TABLE I.- INDEX TO FIGURES 


Figure 

Mode 

Range of Y 

Range of Q , deg 

2 -14 

X 

0(0. 1)0. 9 

0(10)90 and 15 , 45 , 75 

15-27 

Z 

0.1(0. 1)1.0 

0(10)90 and 15, 45 , 75 

28 - 4 o 

z 

1.1(0. 1)5-0 

0(10)90 and 15 , 45 , 75 

41-53 

0 

0(0.1)1.0 

0(10)90 and 15, 45 , 75 

54-66 

0 

1.1(0. 1)5.0 

0(10)90 and 15 , 45 , 75 

67-76 

X 

0(0. 1)0. 9 

0(10)90 and 15, 45 , 75 

77-86 

z 

0.1(0. 1)1.0 

0(10)90 and 15 , 45 , 75 

87-94 

z 

1.5(0. 5)5.0 

0(10)90 and 15, 45 , 75 

95-105 

0 

0(0. 1)1.0 

0(10)90 and 15 , 45 , 75 

106-113 

0 

1.5(0. 5)5 

0(10)90 and 15, 45 , 75 



ooooooon 


TABLE II.- COMPUTATION OF GROUP AND PHASE REFRACTIVE INDEX 


1 

2 
3 

WAIN PROGRAM *** 4 

3 

-I 15 ! 1 '!* INPUT DATA 6 

** (1) NTHtTA = NO. OF THETA USE FOR COMPUTATION (13) 7 

** (2) THETA = ANGLE THETA IN DEGREES (F7.2) b 

(NTHETA OF THEM, ONE THETA PER CARD) 9 

1U 

DIMENSION X( 2UO) , SO ( 2UU ) , SX< 2 0U) , RMU OS (■) ( 2 OU ) , Ri-iUXSO < 200 ) ,R ( 200 > , 11 

* RMUOP ( 200 ) jRMUXP ( 200 ) 12 

DIMENSION TR ( 100 ) , XX ( 100 ) ,TSO< 100) ,TSX( 100 ) , THUOSO ( 100 ) , 13 

* T (iU XSO ( 100 ) » THUO P ( 100 ) * Tl-iU X P ( 100 ) 14 

COMHON/CUf-ilNP/C, S, Y, THETA 13 

COMMON /CUMOUT/Nt A , XX , TSO, TSX, TMUOSO , THUXSU ,T MUOP , Tl'iUXP , TR 16 

RE AD ! 5 > 8 ) NTHETA 17 

8 FORMAT (13) 16 

DO 300 1=1, NTHETA IV 

RE AO ( 5 , 9 ) THETA 20 

9 FORMAT ( F 7. 2 ) 21 

THETAR = THETA*1.745329E-2 22 

C = CUS(THETAR) 23 

S = SIN (THETAR) 24 

C**** COMPUTE V AND A = ( 1-Y**2 ) / ( 1-YL**2 ) , V = 0(0.115 28 

0**4* N = TOTAL NO. OF POINTS AROUND X = 1-Y , 1 , A , AND 1 + Y FOR EACH Y 26 

DU 25u J = 1 ,51 2 I 

N =0 2 b 

TJ = J 29 

Y = <TJ-1.0)*0.1 30 

A = ( 1.0-Y**2) /< 1.0- ( Y#C)**2) 31 

C**** CONVERT A TU 2 DECIMAL PLACES 32 

NA = A* 100 .0 33 

TA = NA 34 

A = Ta*0.01 35 

1F(Y.GT.0.9) GO TO 14 36 

X 1= 1 . 0— Y 37 

X 2 = 1 .U 3b 

X3 = l.U+Y 39 

GU TU 15 40 

14 XI = 1.0 41 

X2 = 1.0+Y 42 

15 IF (Y.GT. 1.0) GO TO 16 43 

C**** NX = NU. UF X VALUES USE FOR CUM PU TAT ION , <X=0(U.1)2 + Y IF Y LESS 44 

C THAN OR EOUAL TO 1, X=0(0.1)10 IF Y GREATER THAN 1) 45 
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NX = (2.U + Y)*10.O + 1,0 46 

GO TO 17 47 

16 NX = 101 48 

17 WRITE (6*18) 49 

18 FORMAT ( 1H1,40X,47HCGMPUTAT ION OF GROUP AND PHASE REFRACTIVE INDEX) 50 

WkITE (6*19) THETA, Y 51 

19 FORMAT ( 1H0,9X,6HTHETA=,F4. 1 , 2X , 2HY= , F3 . 1 ) 52 

WRITE (6,21) 53 

21 FORMAT ( 1H0, 11X* 1HX, 12X, 1 HR* 11X, 5 HMU OS Q, 1 IX, 2 HSO, 1 2X , 7HMU0 PR I M * 7X , 54 

*5HMUXSQ , 11X*2HSX, 12X, 7HMUXPR IM// ) 55 

C**** COMPUTE REFRACTIVE INDEXES FOR VARIOUS X 56 

DU 100 l< = 1 * N X 57 

IK = K 58 

GX = ( T K - 1 ♦ 0 ) J -' 5 0 ♦ 1 59 

GK = SQRT( (Y*S*S)**2 + <2.0*C#(1.0- GX ))**2) 60 

CALL URD(GX,Y,GK ,GSO, GMUOSQ , GMUOP ) 61 

CALL E XORD ( GX , Y , GR * GSX , GMU XSQ , GMUXP ) 62 

X ( K ) = GX 63 

R ( K ) = GR 64 

S 0 { K } = GSO 65 

S X ( K ) = GSX 66 

RMUOSQ(K) = GMUOSU 67 

RMUXSQ(K) = GMUXSU 68 

RMUOP(K) = GMUOP 69 

KMUXP(K) = GMUXP 70 

100 CONTINUE 71 

C**** COMPUTE REFRACTIVE INDEXES AROUND X= 1 , 1 -Y , 1 + Y , ( 1-Y**2 ) / ( 1-YL**2 ) 72 

CALL IN TER‘M (XI) 73 

CALL I NT ERM ( X2 ) 74 

CALL INTERA(A) 75 

IF (Y.GT.U.9) GO TO 20 76 

CALL IN T ERM ( X3 ) 77 

**** REARRANGE SET XX(K) AND THEIR ASSUC1ATE VARIABLES INTO ASCENDING 78 

**** ORDER 79 

20 NM1 = N-l 80 

DO 30 K = 1 , NM 1 81 

KP1 =K+ 1 82 

DO 30 L = K P 1 , N 83 

IF ( XX ( K)-XX(L) ) 30,30,25 84 

25 T E M P 1 = XX { K ) 85 

T E M P 2 = TSU(K) 86 

TEM P3= TSX(K) 8/ 

T L M P 4 = TMUOSU(K) 88 

T E M P 5 = TMUXSO(K) 89 

T t MP6= TMUUP(i<) 90 
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T E M P 7 = TMUXP(K) 91 

TEMPB= TK(K) 92 

TR ( K ) =TK ( L ) 93 

X X { K ) - XX (L) 94 

T S U ( i\ ) = T30 ( L ) 93 

TSX(K) = T S X ( L ) 96 

TMUUSW ( K ) = TMUOSU ( L ) 97 

TMUXSU(K)= TEiUXSU ( L ) 9b 

Tt-iUOPlK) = TMUUP(L) 99 

T».UXP(K) = TEiUXP(L) 100 

XX ( L ) = TEMPI 101 

TSU ( L ) = TEE1P2 102 

TSX ( L ) = TEMP3 103 

TMUUSO ( L ) = TEHP4 104 

TMUXSW(L)= TEEiP3 1 03 

TMUOP { L ) = TEE.P6 106 

TMUXP(L) = TEE1P7 107 

T R ( L ) = TEE'iP B 10B 

30 CONTINUE 109 

h =i no 

OU 200 K = 1 » N X 111 

32 IF ( Ei * G T • N ) 00 TU 44 112 

I F { XX ( Ei ) • OT • X ( K ) ) GO TO 44 113 

IF (X(K)-XX(EU.GT. 0.009) GO TO 33 114 

E,= H+l 113 

GO TU 44 116 

33 IF( AbS(XX(M)-XX(M+l) ) .GUO. 009) GO TO 40 117 

N — Ei+ 1 1 1 B 

GU TU 35 119 

40 WK I TE ( 6 ,43 ) XX ( Ei ) , TK( M ) , TEiUOSO ( Ei ) , TSO (M ) , TEiUOP { Ei ) , TEiUXSO ( El ) , 12 0 

* T S X ( EW , T Ei U X P ( El ) 121 

El=Ei+ 1 122 

GU TU 32 123 

44 WRITE (6,45) X(K),R(K) , RMUOS U ( K ) , SO ( K ) ,REiUOP{K) y RMUXSO(K) , 124 

* S X ( K ) , KEIU XP ( K ) * 125 

43 FUKEiaT ( 1H ,F14.2,4X t E12.5 y F12.5,4X,E12.5,F16.5fF12.5,4X y E12.5, 126 

*F16.5) 127 

200 CONTINUE 12B 

250 CONTINUE 129 

300 CONTINUE 150 

STUP 131 

END 132 

133 

134 

135 
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SUBROUTINE TO COMPUTE REFRACTIVE INDEX FOR X AROUND 1-Y,1,1+Y 136 

137 

SUBROUTINE INTERM (XI) 138 

DIMENSION TR(IOO) ,XX( 100 ) ,TSO< 1 00 ) , TS X ( 1 00 ) ,TMUOSU< 100 ) , 139 

* TMUXSU ( lOO)tTMUOP(lOO) f TMUXP(100) 140 

CUMMUN/CCJMINP/C, S, Y, THETA 141 

COMMON/CUMOUT/N, A f XX,TSOt TSX,THUOSO,TMUXSg,TMUOP»TMUXP f TK 142 

IF(Xl.LE.O.O) RETURN 143 

C*** X VARIES IN STEP OF 0.01 144 

DO 10 L=2,10 143 

T L = L 146 

X={ TL-1 .0 ) *0.01 + (Xl-0.1) 147 

IF (X.LT.0.0 ) RETURN 148 

R= SQRT( (Y*S*S)**2 + ( 2 . 0*C* ( 1 ,0-X ) ) ) 149 

CALL ORD(X f Y,R, SOtRMUOSOtRMUOP) 150 

CALL EXORD(X,Y,R,SX,RhUXSQ,RMUXP) 151 

N =N+1 152 

TR ( N ) = R 153 

XX ( N ) = X 154 

T SO ( N ) = SO 155 

T $ X ( N ) = SX 156 

TMUOSQ(N) = RMUOSQ 157 

TMUXSO(N) = RMUXSO 156 

TMUOP(N) = RMUOP 159 

TMUXP(N) = RMUXP 160 

10 CONTINUE 161 

RETURN 162 

END 163 

164 

165 

166 

C SUBROUTINE TO COMPUTE REFRACTIVE INDEX FOR X AROUND < 1-Y**2 ) / ( 1-YL**2 167 

C 168 

SUBROUTINE I NT ERA ( XI ) 169 

DIMENSION TR ( 100 ) » XX ( 100 ) » TSOI 100 ) , TS X U 00 ) , TMUOSQ ( 100) » 170 

- TMUXSU (100) ,ThUOP( 100) , TMUXP ( 100) 171 

COMMUN /COMI IMP /C , S , Y , THETA 172 

COMMON/CUMOUT/N, A, XX fTSOtTSXt TMUOSQ, TMUXSU t TMUOPt TMUXP, TR 173 

IF (XI. LT. 0.0) RETURN 174 

C*** X VARIES IN STEP OF 0.01 i 75 

DO 10 L = 2 > 10 176 

TL=L 177 

X=(TL-1.0)*0.01 + XI 178 

K = SORT { (Y*S*S)**2 + ( 2. 0*C* ( 1 .0-X ) ) **2 ) 179 

CALL 0 RD ( X , Y , R , SO , RMUOSO , RMUO P ) 180 


li 



CALL EXQRD(X,Y,R,SX,RMUXSQ,RMUXP) 181 

N =iM + l 182 

TR(lM) = K 1 83 

XX { N ) = X 184 

T S 0 ( N ) = SO 185 

TSX ( N ) = SX 186 

TMUOSQ(N) = RMUOSQ 187 

TMUXSQlN) = RMUXSQ 188 

TMUGP(N) = RMUOP 189 

TMUXP(N) = RMUXP 190 

10 CONTINUE 191 

RETURN 192 

END 193 

194 

195 

196 

C SUBROUTINE TO COMPUTE THE ORDINARY WAVE GROUP AND PHASE REFKmCTIVE 197 

C INDEX (MUOPRI AND MUO) 19b 

C 199 

C X = (FN/FJSU 200 

C Y = FH/F 201 

C THETA = ANGLE THETA IN DEGREES 202 

C MUO = PHASE REFRACTIVE INDEX 203 

C MUOPRI = GROUP REFRACTIVE INDEX 204 

C 205 

SUBROUTINE ORD ( X, Y, R t SO , MUOSQ tMUOPR I ) 206 

REAL MUO 1, MU02» MUOSQ, MUO, MUOPP 1 , MU0PP2 , MUOPRI 207 

CUMMON/COMINP/C»S» Y, THETA 2 O 8 

IF(THETA. EQ. 0. 0. AND. X.GE. 0.99999. AND. X.LE. 1.00001) GO TO 39 209 

S01=Y*(1.0+C*C*( 1.0-2.0*X) )+R 210 

S02=Y*$**2 + R 211 

SO = S01/SO2 212 

MUO 1 = Y*( 1 .U+C*#2 ) + R 213 

MU02= SOI 214 

MUOSO = <1*0-X>*{MU01/MUCJ2> 215 

IF (MUOSQ.LE.O.O.OR.ABS(MUOSQ) .LT.l.OE-5) GO TO 40 216 

MUO = SORT { MUOSO ) 217 

MUUPPla < l,0-X)-( (Y*S*S)*( l.O+X) )/R 218 

M U U P P 2 = ( Y * C * C * X ) / ( S 0 * * 2 * ( Y * S * * 2 + R ) ) 219 

MUOPRI = ( 1.0-MUOPP2*MUQPP1 J/MUO 220 

RETURN 221 

39 SO = 1.0 222 

40 MUOSO = 0.0 223 

MUOPRI =10000.0 224 

RETURN 225 
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227 

228 
229 

C SUBROUTINE TO COMPUTE THE EXTRAORDINARY WAVE GROUP AND PHASE 230 

C REFRACTIVE INDEX(MUXPRI AND MUX) 231 

C 232 

C X = ( FlM/F ) SO 233 

C Y = FH/F 234 

C THETA = ANGLE THETA IN DEGREES 235 

C MUX = PHASE REFRACTIVE INDEX 236 

C i*i IJ X P R I = GROUP REFRACTIVE INDEX 237 

C 238 

SUBROUTINE E XORD ( X , Y , R f SX , MU XSO ,MUXPRI ) 239 

REAL MUX1 ,MUX2 , MUX3 , f*1U X4 T MUX24 , MUXSG t MUX ,HUX PP 1 , MUXPP2 , HUXPRI 240 

COMMON/ CUM I NP /C, S, Y, THETA 241 

IF(THETA.EO.O.O) GO TO 60 242 

IF ( ABS< l.O-YJ.GE.O.Ol) GO TO 10 243 

IF(X.LT.O.Ol) GO TO 15 244 

10 I F ( ABS ( 1.0-X).GT. 0.001) GO TO 20 245 

IF (Y.EQ.0.0 ) GO TO 75 246 

SX = 10000.0 247 

MUX SO =1.0 248 

MUX PR I = 1.0 + 1.0/((Y*S)**2) 249 

RETURN 250 

20 SX1 = ( 1.0-X)+Y**2*(X*C**2-1.0) 251 

SX 2 = 2.0*(1.0-X)-( Y*S)**2+Y*R 252 

SX = 2.0*1 SX1/SX2) 253 

25 M U X 1 = 1.0+Y-X 254 

MUX 2= SX1 255 

M U X 3 = U.0-X)*SX2 256 

iiUX4= 2 .0* ( 1.0-X)**2- { Y*S )*#2+Y*R 257 

i-i U X 2 4 = MUX 2*MU X4 258 

MUX SU =(1.0-X-Y)*( (MUX1*HUX3)/MUX24) ^59 

IF (MUXSO.LE.O.O.OR.ABS(MUXSU) .LT.1.0E-5) GO TO 50 260 

MUX = SQKT(MUXSQ) 261 

MUXPP1 = (X*( 1.0-SX) )/(2.0*$X**2) 262 

MUX PP 2 = 1.0+{(Y*S**2)/R)*((1.0+X)/(1.0-X)) ^63 

MUXPKI = ( l.Q+MUXPPl*MUXPP2)/MUX 264 

RETURN 265 

60 IF (ABS ( 1.0-Y) .GE.0.01) GO TO 65 266 

IF(X.LT.l.O) GO TO 49 267 

IF (ABS ( l.U-X) .LT.0.U01 ) GO TO 15 26b 

S X = 2 . 0 269 

MUXSO = 1.0 - X/( 1.0+Y) 270 
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IF ( MUXSQ. LE.U .0.GR.A6S ( MUXSQ ) . LT . I .OE-5 J GO TO 50 271 

MUXPKI = tl.U - ( X*Y) /( 2.0*( 1.0+Y)**2 ) J/SQRT (MUXSQ) 272 

RETURN 2 

65 IF(Y.LT.l.U) GO TO 70 274 

IF( aBS( l.U-XJ.GT. 0,001) GO TO 20 275 

SX = 10000. 276 

MUXSQ =1.0 277 

MUXHkI = 10000.0 278 

RETURN 27V 

70 IF(AbS(1.0-X).GT. 0.001) 00 TO 20 280 

IF(Y.EQ.U.O) GO TO 75 281 

GO T 0 15 282 

75 SX = 1.0 2 85 

GO Tu 50 2 84 

15 SX =0.0 2b5 

MUXS 0 = 10000.0 <i66 

MUX F K I = 10000. 0 267 

RETURN 286 

49 SX =0.0 289 

50 MUXSQ = 0.0 290 

MUX HR I = 10000.0 291 

RETURN 492 

END 253 
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Figure 3*- Variation of \±' vs. X; Y =0-0.9; 0 = 10 
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Figure 7*- Variation of (jl t vs . X; Y = 0 - 0.9; 0 = ^-0° . 
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Figure 8.- Variation of |i' vs. X; Y =0-0-9; 9 = b 5 • 


22 


IjUXfjJR I|1 



Figure 9 *“ Variation of ja* vs.X; Y = 0 - 0.9; 9 = 50°. 
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Figure 11.- Variation of |i f vs . X ; Y = 0 - 0 .9 ; 0 - 70 . 
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Figure 13*- Variation of |i' vs . X; Y = 0 - 0 -9; 9 = 80°. 


27 


tlUXIJRIJI 



28 



^UXf^RIJI 



Figure 15 •- Variation of jjt' vs. X; Y = 0.1 - 1.0; 9 = 0°# 
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Figure 17 •- Variation of |i' vs. X; Y =0.1 - 1.0; 9 = 15 • 
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Figure 19.- Variation of |i' vs . X; Y = 0.1 - 1.0; 9 - 30 . 
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Figure 21.- Variation of |i* vs. X; Y - 0.1 - 1*0 j 9-^5 • 
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Figure 25-- Variation of n* vs . X; Y = 0.1 - 1.0; 9 = 75° • 
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Figure 26.- Variation o: 
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X t : 0.1 I0.1U.0 


H' vs. X; y = 0.1 - 1.0; 6 = 80° . 




^XPRIJI 



Figure 27-- Variation of n' vs . X; Y = 0.1 - 1.0; 9 = 90°. 
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Figure 28.- Variation of |a' vs . X; Y = 1.1 - 5.0; 9 = 0 ° • 


42 


Z TRACE THETA =10 DEGREES 



0.50 1.00 1.50 t.00 £.30 3.00 3.30 4.00 4.30 $.00 S.SO 4.00 

X V = 1.1 ( 0.1 ) 3.0 


Figure 29 .- Variation of fi* vs . X; Y = 1.1 - 5*0; 0 = 10 
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Figure 31 •- Variation of |i' vs. X; I = 1.1 - 5.0; 9 = 20°. 
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Figure 33 •- Variation of |i' vs. X; Y = 1.1 - 5*0; 9 = 40 . 
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Figure 3^*~ Variation of | a T vs . X; Y = 1 .1 - 5 *0 ; 0-^5 
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Figure 36 .- Variation of (i’ vs. X; Y = 1.1 - 5*0; 9 = 60°. 
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Figure 39*” Variation of |i’ vs. X; Y = 1.1 - 5*0 > 0 = 8o° . 
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Figure 40.- Variation of jj.’ vs. X; Y = 1.1 - 5* 0 > 0 - 90 » 
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V = 0.0 ( 0.1 ) 1.0 


Figure 4l.- Variation of p* vs. X; Y 


1 . 0 ; 9 = 0 ^ 
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Figure 42.- Variation of \±' vs . X; Y =0-1.0; 0 = 10 « 
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Figure 44 .- Variation of* ijl* vs . X; Y = 0 - 1 .0; Q = 20°. 
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Figure V?.- Variation of |_l* vs . X; Y = 0 - 1.0; d = 30°. 
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Figure 47*-. Variation of |i’ vs. X; Y = 0 - 1 . 0 ; 9 = ^5°* 
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Figure 48 .- Variation of n' vs. X; Y =0-1.0; 9 = 50°. 



62 


. WRIJ1 



Figure 49.- Variation of (i 1 vs . X; Y =0-1.0; 9 = 60° . 
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Figure 51 •- Variation of [i 1 vs. X; Y =0-1.0; 9 = 75 ° • 
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Figure 5 2 Variation of |i' vs . X; Y = 0 - 1.0; 0 - 80° . 
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Figure 53-- Variation of n’ 


vs . 


X; Y 


0 


i.o; e = 90°. 
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Figure 55*- Variation of i± l vs • X; Y = 1 .1 - 5-0; 0 = 10 • 
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Figure 57*- Variation of ji 1 vs . X; Y = 1.1 - 5*0; 0 = 20 . 
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Figure 58 Variation of n' vs . X; Y = 1.1 - 5»0; 0 = 30 
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Figure 59*“ Variation of p.’ vs. X; Y = 1.1 - 5*0} 9 = ^0° 
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Figure 60.- Variation of n' vs . X; Y = 1.1 - 5-0; 0 = ^5° • 
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Figure 6l.- Variation of |f vs. X; Y =1.1 - 5.0; Q = 50°. 
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Figure 63-- Variation of ju' vs. X; Y = 1.1 - 5»0; 9 = 70°. 
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Figure 64.- Variation of |a* vs. X; Y = 1.1 - 5*0; 9 - 75° • 
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Figure 65 .- Variation of |i' vs . X; Y = 1.1 - 5*0 ; 0 = 80°. 
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Figure 67-- Variation of |i 1 vs. X; y =0; Q = 0° - 90°. 
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Figure 69 .- Variation of (i 1 vs . X; Y =0.2; 0 = 0° - 90° • 
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Figure 72.- Variation of p.’ vs. X; Y = 0.5.5 0 =0° - 90° • 
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Figure 7^«- Variation of |F vs. X; Y = 0.7j 0 = 0° - 90° • 
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Figure 75 •“ Variation of |i' vs. X; Y = 0.8; 9 = 0° - 90° • 
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Figure 77*- Variation of p' vs. X; Y = 0.1; 0 = 0° - 90 ° • 
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Figure 8l.- Variation of p' vs. X; Y = 0.5; 0 = 0° - 90° • 
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Figure 82.- Variation of |i' vs. X; Y = 0.6; 9 - 0° - 90° • 
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Figure 93*” Variation of |i* vs. X; Y - 0 - 0° - 90°. 
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Figure 95*“ Variation of n' vs . X; Y = 0; 0 = 0 C 
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Figure 9 6.- Variation of |i' vs. X; Y =0.1; 9 =0° - 90° • 
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Figure 99*- Variation of [x' vs . X; Y = 0.4; 0 = 0° - 90°. 
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Figure 100.- Variation of |F vs. X; Y = 0*5; 6 = 0 ° - 90°- 
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Figure 103 •- Variation of n T vs. Xj Y = 0.8; 9 = 0 ° - 90° • 
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Figure 10^.- Variation of |i' vs . X; Y = 0.9; 0 = 0° - 90°. 
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Figure 107 •- Variation of |i' vs. X; Y - 2.0; 6 - 0° - 90 • 
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Figure 111.- Variation of |a' vs . X; Y = ^-.0; 9 = 0 ° - 90°. 
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Figure 112.- Variation of n’ vs . X; Y = 4.5; 9 = 0° - 90° • 
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Figure 113.- Variation of n’ vs. X; Y = 5 . 0 ; 6 = 0° - 90°. 
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